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THE DETERMINATION OF THE MASS OF A PLANET FROM THE 
RELATIVE POSITION OF TWO SATELLITES. 


By Pror. AsapH HALL, Washington, D. C. 


OnE of the easiest methods of determining the mass of a planet is from the 
observations of its satellites. The two principal elements that are needed are 


the periodic time of the satellite about its primary, and the semi-major axis of 
its orbit. Thus, if M7 and m are the masses of the sun and the planet, and +, 
and z are the periodic times of the planet and the satellite, we shall have very 
nearly, 

a® 

a>’ 


-2 
. 


a, and a being the mean distances of the planet from the sun and of the satellite 
from its primary. The above is the simplest form of this equation, and should it 
be necessary in any case to consider the elliptical figure of the planet, and its mass 
when added to that of the sun, these small corrections can be applied easily. Of 
the two principal elements the periodic time of a satellite can be found with suf- 
ficient accuracy by comparing observations of it made at distant epochs, and 
hence the error arising from this source may be considered as vanishing. But 
with respect to the semi-major axis the case is different. For the determination 
of this element we have to depend on micrometrical measurements, and several 
sources of error present themselves. Thus the value in arc of one turn of the 
screw of the micrometer must be known very exactly, since this value will fre- 
quently be multiplied by a large factor, and of course an error in the value of 
the screw will also be multiplied by the same factor. With care on the part of 
the astronomer this error may be reduced to a very small quantity; and yet it is 
surprising how often such errors remain in micrometers that have been a long 
time in use. But there is another error in micrometrical measurements that is 
not easily got rid of, and this is the personal error of the observer. This error 
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will sometimes amount to half a second of arc, and since it is not eliminated by 
the repetition of observations, each observer will have his own value of the mean 
distance of a satellite, and hence his own mass of a planet. This personal error 
may be still greater, if, as is generally the case, we have to do with the disk of 
a planet. In this case we must bisect the disk with one of the wires, or else 
measure from the limb of the planet, and both operations are difficult and tend 
to increase the constant error of the measurement. 

On account of this difficulty Struve, the present Director of the Pulkowa 
Observatory, has proposed that observations be made of the relative positions of 
two satellites for the purpose of correcting their elements, and in this way de- 
ducing the mass of the planet. Mr. Backlund, one of the astronomers of the 
Pulkowa Observatory, has already begun a series of measurements of the rela- 
tive positions of the satellites of Jupiter with a heliometer, and it is said that 
these observations are two or three times as accurate as those made by reference 
to the limb of the planet, and moreover that the complete observation of the 

relative position of two satellites requires only a 
third part of the time that is necessary to de- 
termine the position of one of them with respect 
to the center of the planet. Mr. Backlund has 
published an outline of this method, which is 
theoretically very simple. “Let P be the place 
of the planet projected on the heavens, // the 
pole of the heavens, S, and S, the satellites, and 

aa denote the angles and sides of the spherical tri- 
angle, as in the figure. 

The angles of position and distances are designated generally by f and s, 
with proper subscripts, and the relative angle of position and the distance of the 
two satellites by // and a. Then 


angle S, PS, =p. —/1; angle S,S, —7’, 


and the spherical triangle S, PS, gives the equations, 

sin sin —7z’) = sin s, sin —/P,) 

sin cos —2z’) = sins, cos s, —cos $s, sin cos (Py —/P}). 
These are exact formulas, but as s, and s, are always less than 12’ in any real 
case that we have, we may generally put 

J, = ond 

COS = COS =I, and sin SiN S,=Sq, sino=—<a, 

and hence we get 
asin (~,; —Z) = Sq sin (py 


a cos —Z) = —S_q cos (Py 


3 
ace. 
4 
Ay 


HALL. DETERMINATION OF THE MASS OF A PLANET. 


and from these we have 
=z = s, sin —sy sin py 
Z = $1 COS Py —Sy COSPy. 

Now s, sin f, and s, cos f, are the rectangular co-ordinates of the first satel- 
lite with respect to the center of the planet, and the preceding equations become 
where « and y are the relative co-ordinates of the satellites; a result that might 
have been anticipated, since the above assumptions reduce the work to that of 
plane co-ordinates. Should it be necessary in any case we may use the rigor- 
ous formulas; or, what is easier, we may expand them and retain terms of the 

lower orders. In the present case we have therefore 

dx = dx, ax, 

dy = ay ay. 
The differentials on the right hand side can be expressed in terms of the differ- 
tials of the elements of the orbits of the two satellites, and of the differential co- 
efficients with respect to their elements; and since the left hand side will be 
known by comparing the computed relative position with the observed, we have 
the means of forming the equations of condition for correcting the elements of 
these orbits, and therefore for determining the mass of the planet. Thus dv, 
and dx, will be of the form 


1 


where a,, ¢,, P;, @, 1,7; are the elements of the orbit of the first satellite, and 
Py, Ny, Jo are those of the second satellite. The equations of condi- 
tion therefore that result from this method will each contain twelve unknown 
‘quantities ; and it will be necessary to arrange the observations so that there may 
be no indetermination in the solution. For this purpose the satellites should be 
observed when they are on different sides of the planet as well as when they are 
on the same side, and in a variety of positions. 

In the common method of referring a satellite directly to the center of the 
planet we have in the equations of condition six unknown quantities. Now the 
labor of solving a set of equations of condition increases rapidly with the number 
of unknown quantities. If 7 be this number, then the number of coefficients in 
the normal equations is found by the summation of an arithmetical series to be 
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The number of auxiliary quantities that have to be formed in the successive 
elimination of the unknown quantities from the normals is given in the following 
table, with their first and second differences. 


No. 2 
22+ 2 
2 
—t+o 
2 
+1 
2 
727 —37+ 0 
> 4 - 
27—5i-+- 4 
2 
3 
12 —7i-+- 10 
2 
i+ 4 
22? —ot +18 
2 
etc. 


The law of this series is evident. The sum of 7—1 terms of the series is 


(@—1) (¢+2) 22 (¢—1) ¢(@—1) 
2 2 2.3 
(¢—1) (¢+ 4) 
2.3 
Hence the whole number of quantities to be formed in this solution is 
(¢—1)¢(¢-+4) , _ (¢+1) (#45) 
2.3 2 2.3 
In the case of 76, we shall have 77 quantities to form, and when 7=— 12, we 
shall have 442 quantities to form. The labor of solution is nearly six times as 
great in the new method, but this is a consideration that is not perhaps entitled 
to much weight. Again, when the number of unknown quantities is increased, 
the probable error of a single equation is generally increased, but this increase 
will be small when the number of observations is large; and if the observations 
be well arranged, this question may be left out of consideration. On the other 
hand, the gain in accuracy seems to be substantial, and the new method prom- 
ises a real advance, especially in the case of. the satellites of Jupiter and Saturn. 
Another advantage is that it is a different method of making the observations, 
and in this way it offers an independent check on the old methods, and may be 
the means of eliminating some constant errors. 
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ON THE LUNAR INEQUALITIES PRODUCED BY THE MOTION OF THE 
ECLIPTIC. 
By Dr. G. W. HiLt, Washington, D. C. 


Tuis subject has been treated by Hansen,* and more recently by Sir G. B. 
Airy and Prof. J. C. Adams.+ Hansen’s discussion is accommodated to the 
peculiar system of co-ordinates he employs, and the two later writers do not 
consider the inequalities in longitude. Hence an investigation, giving the ine- 
qualities of the latitude and longitude, at first, in the literal form, may be of 
value. The procedures employed are very similar to those of Pontécoulant, and 
are doubtless not as direct as might be imagined. The paper was written as long 
ago as 1867. 

I. 

Expressed in the ordinary notation, when the co-ordinates are referred to 

fixed planes, the lifferential equations of motion are 


at? ax’ 
dy’ 
a@Z aR 


Since the directions of the axes are arbitrary, let the axis of .Y be directed 
toward the ascending node of the moving ecliptic on the ecliptic of 1850; and 
let the axis of Z be perpendicular to the latter plane. Taking now another sys- 
tem of co-ordinates, x, 7, and s, such that the axis of + has the same direction 
as that of Y, but the axis of < is perpendicular to the moving ecliptic, let 
xz (t— 1850) be the inclination of the moving ecliptic to that of 1850; then, 
neglecting quantities of the order of =”, these equations exist, 

A = 
Y— y —z(¢— 1850) z, 
Z=2-+ (t¢— 1850) 7. 
The differential equations of motion, expressed in terms-of the second system 


of co-ordinates, are 
aR 

at? rs ax 
aR 
aR _ day 
at? rs" dg dt 

* Darlegung etc., Arts. 175-178. 

t Monthly Notices, Vol. XLI., pp. 264, 375, and 385. 
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Denoting the true longitude of the moon by 4, from these may be derived the 
two 


— 
at ah 


In this discussion all terms involving the solar eccentricity and parallax will be 
neglected. Let £ denote the moon’s mean angular distance from a point go” 
behind the ascending node of the moving ecliptic on the ecliptic of 1850, or 
£- ¢-- ut—TIl— go°. For simplicity, the semi-axis major of the lunar orbit 
will be made equal to unity, and, as usual in the lunar theory, 7 will be written 


for —. Also let ¢ and z denote respectively the true and mean angular distance 
n 


of the moon from the sun. 
With these restrictions and notation 


aR "AR 
2 -4R+- 2m - a”, 
9 
m? 
R 7 s*) cos 2¢ + r?— 327], 
aR ; 
adh 2 


If the symbol ¢ prefixed to any quantity denote that part of it, in its develop- 
ment in series, which is multiplied by the first power of z, the equations for de- 


termining 07, 0/4, and ds are 


» ‘ » 
yor dR vas — sdy 
=5 == 40R 2m 0. + 2-———-, 
+25, (rar — 282): 0. + 2 = 
a ‘a e ash %. 
ors m2 | —2-—. 


In these equations terms multiplied by the square and higher powers of the in- 
clination of the moon's orbit are neglected; and, since dr and d4 are multiplied 
by the first power of this quantity, this involves the neglect of terms such as sdr 
and 204. ‘For the same reason all higher powers of < than the second have been 
omitted in AX. 

These equations suffice to determine the inequalities we seek; but for a 
term of long period in d¢ it will be more commodious to employ another equa- 
tion. We have 


a?r », ar? (dz sdr\?. yds — sdy 


dz? av? lad r dz ) dr ad” 
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or 
( 1dz sdr\? R z yds — sdy. 


Taking the variation with respect to z, and then multiplying by 7?, 


» dh ah, ra? or —- d?r.dr | aR aR, 
| 2 >, | — 2—. dA— 2— dz 
| ae ae ar? | | ah az | 
| z ar Os 3yuror | zyds—sdy | 
de) ae rae n ade 
But 
ad? ror uror dR x vydz—sdy 
3 = 120R - 6m { 0. - ad” +6-- 


subtracting this 


= = = 502 
de ae dz? ar? dy az ade raz 
aR , aR zx yds — sdy 
sos —70R - Or — 3772 0. ar 


In determining dz we shall stop at terms of the order m —, and shall neglect 
nN 


all terms multiplied by powers of the lunar eccentricity ¢ higher than the first. 
In 04 we shall neglect ¢ altogether; and, since the inequalities in the lunar par- 
allax resulting from dy are insensible, dv will be determined only so far as it is 
necessary to the determination of 04. Let ¢ denote the moon’s mean anomaly, 
and 7 its mean argument of latitude, or 7 =¢-xt— &. In applying the last 
equation to determining the co-efficient of sin (f — 7) in dA to terms of the order 


of 7 —, (where 7 denotes the same function of the inclination as it does in Ponte- 
n 
coulant’s 7héorte Analytique), it will be necessary to compute each member to 


terms of the order of m?;—. But 7dr is of the order of 7-, consequently 
a?ror 
yr? in the term which has £—y for its argument, is of the order of m47—, 


ar . ‘ 
and thus may be neglected; moreover qe is of the order of m?, and hence 
didr.or) . 


pr «CS of the order of m4; “ in the term having the same argument; this 
may then be also omitted. 
With these simplifications the last equation becomes 
ar a(sdz) 
ar ae rae az 


x yds — sdy 
m? (1 + cos 2¢) — 4-—— 


adid.oh adh? adsd.dz 
at ae ag 


> 


SOS 


nN da” 
— 3m? (1 + 3 cos 2¢) ror — 7 —- 0A — 32m é. 
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If, for brevity, we write 


¢ 
A 
B : — 2m? — 6m? cos 2¢, 


C =6m?r? sin 2¢, 
D —3m?*r? cos 2¢, 


E 3m” sin 2¢, 
zdy 
=2- —6m? (1 + cos 2¢) sds + | sin 2¢. sds az, 
ae 
"= + 2 3202 + 3m? f sin 29. 


the term 2m | ErordZ in the equation for dr being omitted as not giving any 
terms which we wish to preserve, and it being sufficient to put 4 — 1, and 


as 


5 1, where the latter multiplies 707 in the equation for 04, the three equations 
¢ > 


become 
+ Ads 
az? 
——,- + ror -+- 2m Doidz 
ar F 
r? 2 ror + [ Dos — Eror| de = U". 
ay 


To the degree of approximation we desire 


I I9_ . 
I m? -- (m? cos 2t, 
6 6 
II, 
A=e-+at - ( — 59 m | sin 2r. 
8 


Also (Pontécoulant, Théorie Analytique, Tom. IV., pp. 216, 226) 


2 5 Ic 13 
+ 942 m m*) cos 2t 
2" 32 16 (3 2 6 
(45 m | e cos (27 — $) - cos (27 + $). 


to 


(3 16 


32 
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From y sin (A— we derive 


2 
m 19 I 
I= [ cos cos (2 — 2r) cos + &), 
2 - - - 
u=—(2-* sin ¢ — sin — 27) — 2e-sin — ¢). 


A, sin + Ay sin — 27) + Ag sin + 27) + A, sin — 47) 
+A;esin (<—€&)+Age sin (+ €) sin — 27 + + Age cos (£ — 27 —§) 

| 

On substituting this expression in the first of the three differential equations, the 

following equations result for determining A,, Ag, etc., 


+ Age sin (€—2t — + A, sin + 27 + +4, sin (€ — 47 + &) 


A,+ ( 3n2+ 12 137 ns | (A,+A;\=—2- 
2 32 16 2 4 12 , 3 
— (8 — 12m) Ag + (3 | A, O, 
— 8A,+ 3 A, O, 
(3+3 ] +45 wes 
— (3 —2 m3) (3+3m] A, - 2, 
4 A, + 3A. + 45 4 957 

5 A, 5 [45 6a m? | A, O, 
— (3 — 8m) A, +3m? A, +34, + (45 m ue | A, fe) 
— (3 — 8m) A, +3 m? A; +34, - a m? A, O, 

3 14 ‘ 
— 1§ A, + + A, O, 
_4 
O. 
By the solution of these, this expression of dz is obtained. 
= — (44-24! 3395 m ) sin 
(3 2 N 


pa : 
{ 
j 
ii 
Let 
id 
4 
i fm 
} 
- 
AWN 
I 25 , 1843 
if 
| 
iy 


10 HILL. ON LUNAR INEQUALITIES. 
I, 121 
— + } —sin + 27) 
4 96 
+ (€—4r) 
32 
16 
+ E m* 109 | e—sin (€—§&) 
4 90 
2 43) 
2 
— —m e —sin 
3 8 n 
3m eZ sin (€ —2r-+6) 


The value of z ( Théorie Analytique, Tom. IV., pp. 237, 244) is 


s=7 + ) + sin (27 —7) 


217 3 a7 + 
3m } sin (27 


whence, by multiplication, is obtained 


303 = — [ 2+. 7 —cos (€ —7) 
( + me} 77,008 €—y—2r 


Iz 
+ -y-cos 27 


[TO BE CONTINUED. ] 


- 
— | 2m-14+-42 | e=sin (£+ 27 —&) 
4 n 
I 9) ~2r—é 
- 
-~m~1— = | e—sin — §) 
4 24 
I a . ~- 
- 
——¢—om + 2t 
- 
+—e-sin (€—47-+$). 
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RECTILINEAR ASYMPTOTES. 
By Pror. E. B. Smitu, Richmond, Va. 


1. THE radius-vector from the origin to a point at infinity is parallel to the 
asymptote containing that point. Hence, if from the equation to any such 
radius, y= yx, we substitute for y in the algebraic equation to a curve, ¢ (1, 7’) 
=o, and from ¢ (x, 4x) =o ascertain a value # corresponding to +=, we 
obtain the inclination of the asymptote to the 2-axis. 

2. The intercept of the asymptote on the y-axis is equal to the difference ,3 
between the ordinate of the asymptote and that of the radius-vector = py, 4, 
for any common abscissa. Or, as = — py 


x x 
since in ¢, y=yex is the ordinate of the curve which becomes the ordinate of 
the asymptote when the abscissa becomes infinite. Thus, the equation to the 
asymptote is y= py Jy. 

3. It is obvious that no terms in ¢ under the highest or x" degree enter 
into the determination of finite values of #. If there be infinite values, that is, 
if there be asymptotes parallel to the y-axis, these can be found by inspection 
of the co-efficient of the highest powers of y in ¢. 

4. If ¢ be rational, since only finite values of ,3,, are possible for a finite 
asymptote not parallel to the y-axis, in the investigation of these, besides the 
terms of the 7" degree in ¢ (x, x) only those of the (7— 1)” degree are in 


when 


G(X, Ur 
volved. Let these terms be respectively and then 
=% and fig, is found from f=0 and Bo =— 77 
dn 


5. If the coéfficients of the terms of the two highest degrees in ¢ be none of 
them zero except that of 2", one value of w, is Oo and a line parallel to the 
4-axis is an asymptote. If, at the same time, /, has y» for a factor, that is, if ¢ 
has no pure x term of the degree x — 1, then ,4,, =o, and the -+-axis is itself 
the asymptote. The conditions for this result are then, the absence only of the 
pure + terms of the two highest degrees in ¢g. If, also, when x is not less than 
4, the two lowest powers of « (namely, the first and the zero power) be 
wanting, the a-axis is also a tangent, and thus a tangent and an asymptote coin- 
cide. If x= 3, the cubic consists of a conic and the «+-axis, and the latter may 
be regarded as the limiting form of a coincident tangent and asymptote. 

6. If g be a homogeneous 7” the origin will be a multiple point of the order 
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n, and the different values of » will each be the same for every value of x, and 
hence the equation is that of a system of x right lines through the origin, each 
of which may be considered the final form of a coinciding tangent and asymptote. 

7. The general condition for such a coincidence may be thus stated: 

If y= 4p is a tangent to g(x, y)=0, at r=, then ¢ y)= 
— ¥ —- $1 + — 4)? Gg=0, and therefore (x—)? ¢,=0, or the 
remainder on division of ¢ by y—yp, x — 3, equated to zero, will yield the 
abscissa of contact as one of a pair of equal roots. 

[It should be observed, with reference to article 4, that, if fo have equal 


roots, 4 will be o when w=. In this case there will be no asymptote at 
au 


a finite distance unless /,; is then alsoo. In that case ,3, will be determined 


9 


from the relation ,4,, =lim - [4= and so on. 


Again, with reference to article 7, it is to be noted that y — nu, «— 3, =O 

is tangent to ¢g at r=; but it will not be also an asymptote unless ¢, is three 


2 


degrees lower than C1, for the limit of the ratio (x — k) 2 ; [x = 90 ], can not 
< ‘ r 
1 


otherwise be zero.—Ep.] 


CONSTRUCTION OF PERSPECTIVE PROJECTIONS. 
By Pror. Wm. M. THornton, University of Virginia. 


THE common experience of draftsmen as to the laboriousness and tedious- 
ness of the ordinary processes for the construction of perspective projections 
has induced geometers to introduce numerous improvements into these pro- 
cesses. The objects aimed at in these improvements have been the reduction of 
the amount of purely constructive drawing, which in the ordinary processes 
largely exceeds that which appears in the net result or picture; the restriction 
of this constructive work to a narrower space, with a view to the enlargement of 
the picture which could be constructed on a given sheet; and the increase of ac- 
curacy in the determination of the outermost points of the figure, which owing 
to the obliquity of the intersections frequently left much to be desired. The 
improvements have been effected by using plans and elevations constructed on 
different sheets from that which was to contain the picture; and by employing 
different scales for the various sets of lines used in the construction, so as to 
dispense with the employment of remote vanishing points. The works of Ad- 
hemar exhibit these improved methods in perhaps their best form. 

The process explained in this note is that which I habitually expound to 
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13 


my students in engineering. After careful comparative tests of the various pro- 
cesses in the drawing-room, I have concluded that it offers the highest guaran- 
tees of accuracy and permits the most rapid and facile execution of the draw- 
ings. As will be seen, it is an attempt to effect an alliance between purely 
graphical methods and arithmetical processes, the draftsman being left free to 
use that method at each point of the work which seems to lead most directly 
and surely to the desired result. 


y P x 


I refer the points of the original to an orthogonal co-ordinate system of 
three planes. The x-plane is the perspective plane; the y-plane is the vertical 
plane at right angles thereto through the center of projection; the s-plane is the 
horizontal plane through the center of projection. The distances of P from 
these planes are x, y, s; the perspective is referred to the axes traced by the 
yand s-planes on the z-plane; and the distances of z, the perspective of /, 
from these axes are 7, %. The distance of the point of view from the perspec- 
tive plane is @. The figures show the relations of these parts [1] in plan, [2] in 
elevation, and furnish at once the fundamental formule 


V 


é 


It is easy always so to adjust the scale that d shall be 100 or 1000. 

These formulae may be used to compute both co-ordinates of z; or to 
compute one co-ordinate, the point being then located by intersection on a per- 
pendicular, a diagonal, or some line of the figure already constructed. 

In particular they give for the vanishing point of a set of parallel lines whose 
common azimuth is « and declivity 0 


4 —d tana, tan 0 seca. 


> 


Such points furnish invaluable tests of the accuracy of the delineation. 
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SOLUTION OF THE QUARTIC EQUATION, 
x4+-Ax+B=0. 
By Mr. A. M. Sawin, Evansville, Wis. 


From Harriot’s laws we easily derive the root relations 
t=—b, x=}[a+6 3a?+ 20d 36?)], 
whence 
«4+ ab?-+ + + a2b?-+ ab*) =0, 
A = + ab? + 68, 
B= ab — a?b?— ab >. 


Put for a—éand for a— we get 


y?— y* = 168. 
Eliminating y we get the cubic resolvent in the form 

26— 4 Bs2?— A? =—0, 
from which < can be found. We then find y from the relation 


52 
and calculating a and 6 we can form at once the four roots. The form of the 
equation is interesting from its occurrence in occasional problems in the resist- 
ance of materials and in hydraulic engineering. 

[The root combination here used was employed first by Bette in an obscure 
dissertation (Neue Auflisung der biquadratischen Gleichungen, 1854). The form 
of resolvent obtained is that due to Descartes. —Eb. | 

APPROXIMATE EXTRACTION OF ROOTS. 
By Mr. Henry Heaton, Atlantic, Ia. 


‘ root of a number is one 


of its x equal factors, the average of 7 nearly equal factors will be an approx- 


Tus method is based on the remark that as, the 


imation to that root. For example, 2 is a first approximation to the fifth root 
of 40, and since 
2x 


(2--27+2+ 23+ 25) 


to 


is a still closer approximation. This may be made the basis of a new approx- 
imation, and so on. More generally, if «-—- be the x root of R, 
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differs from the true root by a small quantity of the second order only. Since we 
can always take for x the nearest whole number to the root considered, / will 
be a proper fraction less than $, and the error committed will be less than 
Sx’ 
that is, less than 3 (%—1) units in the sz place of decimals, # being the num- 
ber of integral places in x. 
The principal value of the method will be found in extending the limits of 
application of ordinary tables of roots. For example, to find 


3.141593 
we take from a table of cubes 1.46 as the first approximation. We find 


3-141593 1.46? 1.47382; 
€= 1.46061 


within 2} units of the last place. The true root is 1.46459. The next approx- 
imation would give nine places, and so on. 


MR. JAMES GLAISHER’S FACTOR TABLES AND THE 
DISTRIBUTION OF PRIMES. 


By Pror. W. W. JoHNnson, Annapolis, Md. 


THE publication of Mr. James Glaisher’s Factor Table for the Sixth Million* 
marks the completion of a work of extraordinary extent, and of the highest 
value in the theory of numbers,.to which France, Germany, and England have 
contributed equal portions of the tireless zeal of the skillful calculator; namely, 
the tables which enable us to resolve into its prime factors every number less 
than 9,000,000. The present volume is noteworthy also as containing, in the 
form of an introduction of more than 100 quarto pages, a very able memoir 
upon the distribution of primes, embodying in tabular form the results of pro- 
found and laborious research. 

The history of the factor tables is briefly as follows: The subject of the 
construction of extensive factor tables was agitated to a considerable extent 
during the latter part of the last century. Euler published a memoir in 1774 
on the mode of construction of a factor table to extend to a million, and four 
years earlier, Lambert, in a work containing a description of a method of form- 
ing such a table, gave a folding sheet on which were given all the simple factors 


* Factor Table for the Sixth Million, containing the least factor of every number not divisible by 2, 3 or § 
between 5,000,000 and 6,000,000. By JAMES GIAISHER, F. R.S. London: Taylor and Francis, 1883. 
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of numbers not divisible by 2, 3 or 5 from 1 to 10,200, and he urgently appeals 
to calculators in the hope that some one will extend this table to 102,000, or ‘‘ if 
he will make himself an immortal name,” to 1,020,000. But it was not until 
1810 that a complete factor table of the first million appeared. This was the 
Cribrum <Arithmeticum of Chernac, which gives all the prime factors _ of 
every number not divisible by 2, 3 or § up to 1,020,000. Chernac speaks of 
this as the labor, not of months, but of years. Three years later, Burckhardt 
published his table for the second million, or more strictly from the limit of 
Chernac to 2,028,000. This gives the least factor only, which obviously effects 
an enormous saving in bulk, while impairing the usefulness of the table only to 
a slight extent, the process being of course to divide by the tabular factor and 
enter the table for the first million with the quotient. To facilitate the divisions, 
a table of the first nine multiples of all primes up to 1423 was given on the first 
page. In 1816 Burckhardt published the continuation of this table to 3,036,000, 
and in 1817 he published the first million in a table uniform with his other vol- 
umes. He states that this table is based upon a comparison of a manuscript by 
Schenmark, in the possession of the library of the Institute, with Chernac, the 
latter being found correct in nearly every case of discrepancy. Burckhardt ex- 
presses his willingness to prepare the manuscript of the fourth, fifth, and sixth 
millions, if the sale of the three millions, which were also published in a single 
volume, shall appear to the publisher sufficiently favorable. 

In 1850 Gauss wrote to Zacharias Dase, who was possessed of a great 
natural talent for computation, to repeat, as he says, in writing what he had 
before told him by word of mouth. He refers in this letter to a manuscript of 
the fourth, fifth and sixth millions which had been presented by Crelle to the 
Berlin Academy, and adds that he doubts not that this will sooner or later be 
published; he therefore urges Dase to undertake the calculation of the next four 
millions. In 1860, through the support of the patrons of science in his native 
town, Hamburg, Dase was enabled to devote himself entirely to carrying out 
Gauss’s project, but in 1861 he died suddenly, leaving the seventh million com- 
plete and the eighth nearly complete, as well as a great part of the work for the 
ninth and tenth millions. The seventh million volume, in the preface of which 
Gauss’s letter is printed, was published in 1862 by the Committee of the Dase- 
Stiftung. Dr. Rosenberg, of Hamburg, undertook the continuation of the work, 
and the eighth million was published in 1863. The ninth million, described on 
the title page as by Zacharias Dase and completed by Dr. H. Rosenberg, fol- 
lowed in 1865; a short preface signed ‘‘Das Comite der Dase-Stiftung”’ states 
that the tenth million is near completion. This million has, however, never been 
published. 

In the report, for 1878, of the Committee on Mathematical Tables to the 
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British Association for the Advancement of Science, it is stated that, the com- 
mittee deeming it very desirable that the gap of three millions between 
Burckhardt’s and Dase’s tables should be filled up, Prof. Cayley, the chairman 
of the committee, communicated early in 1877 with the Berlin Academy, with 
regard to the probability of the publication of the manuscript of Crelle, and 
learned that it had been examined and found to ‘be so inaccurate that ‘‘the 
Academy was convinced that its publication would never be advisable.’’ The 
calculation was then at once commenced by Mr. James Glaisher, a member of 
the committee, with the assistance of two computers. The table for the fourth 
million was published in 1879, with an extensive introduction largely historical, 
from which the above summary is chiefly taken. The fifth million was pub- 
lished in 1880, and the printing of the sixth million was commenced in that 
year, but the elaborate character of the introduction has delayed the appearance 
of this volume until 1883. 

It has been learned that the manuscript of the tenth million, commenced by 
Dase and completed by Rosenberg, was presented in 1878 to the Berlin Academy 
by the widow of Dr. Rosenberg. 

The arrangement of the tables is uniform for the nine millions. Owing to 
the rejection of the multiples of 2, 3 and 5, the last two figures of the numbers 
retained recur in a cycle of 300, there being 80 numbers in 300 retained. The 
quarto page admitting of 80 lines, the terminal figures of the argument are 
printed at the side of the page in three groups, each corresponding to 100 num- 
bers. The columns are thus divided each into three sections, each section 
being headed by the two figures of the argument preceding those found at the 
side of the page. The entries in the columns are the least factors in the case 
of composite numbers and bars in the case of prime numbers. The mode of 
construction is due to Burckhardt, but being very briefly described by him is 
explained in some detail by Mr. Glaisher in the introduction to the fourth million. 
It is evident that, were all the numbers retained in the table, the multiples of 
any prime f would occur in certain diagonal straight lines, and therefore could 
be determined mechanically. The rejection of the multiples of 2, 3 and 5 dis- 
locates this arrangement, but Burckhardt perceived that the arrangement of 
the multiples in the first # columns would still be repeated in the next / 
columns, and so on in each successive block of f columns. He therefore pre- 
pared his manuscript on sheets having the terminal figures printed at the side, 
and containing 77 columns, a multiple of 7 and 11 being chosen, so that the 
entries of the factors 7 and 11 would recur in the same position on each shect. 
Thus the factors 7 and 11 could be printed upon the sheets, there being 880 
entries of 7, and 480 entries of 11 on each sheet. This not only saves a great 
amount of labor, but ensures the accuracy of these 1360 entries on each sheet. 
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In like manner the multiples of 13 recur in the same positions in successive 
blocks of 13 columns. If a piece of paper of the size of 13 columns of the 
sheet, having the spaces which correspond to the multiples of 13 in the first 13 
columns cut out, be laid upon these columns, the entry 13 may be made in all 
of the spaces of the sheet which appear through the holes unless they already 
contain a7 oran If. The screen thus constructed, which is called a ‘‘sieve,”’ is 
then placed upon the next block of 13 columns and the entry 13 made through 
the holes, and so on. A sieve for 17 is next used in a similar manner, and the 
entry 17 made in those spaces which appear through the holes and are not 
already occupied by 7, 11 or 13. Sieves were thus constructed for all prime 
numbers up to 307. A little consideration will show that each sieve will have 
80 holes, one in each of the 8o lines, and this was a useful check in their con- 
struction. In practice Mr. Glaisher found it more convenient to employ for the 
small numbers sheets consisting of several sieves in juxtaposition (thus the 13,s 
were entered bya sieve of 78 columns, containing 6 copies of the sieve), and 
for the large primes the sieves were constructed from several sheets which were 
not fastened together, but numbered and used successively, each sheet of the 
sieve, of course, in general overlapping two of the sheets of the table. 

For primes greater than 307 another method was employed. It is obvious 
that every number for which a prime greater than 307 is the /east factor, and 
within the limits of the table, is the product of two primes. To determine 
these multiples, in the case for example of 311, we have therefore to multiply 
311 by successive prime numbers. The first multiple occurring in the fourth 
million was 9649 » 311 = 3,000,839; the series of primes commencing at 9649, is 
9649, 9061, 9677, . . . , the differences being 12, 16, 2, . . . ; adding to 3,000,839 
successively 12> 311, 16311, 2*311,... Wwe have the multiples required. 
The multiples of the larger primes were thus prepared in books, from which 
the entries upon the sheets were made in the spaces indicated, after the sieves 
had been used. These spaces should, of course, all be found vacant, and this 
afforded a verification of the correct placing of the sieves. 

The introduction to the fourth million contains a list of primes up to 
30,341, with their differences. 

Considered solely with reference to the frequency with which the results 
may be required, it might be thought that a factor table should not rank high 
among mathematical tables for utility, but this measure of comparative value is 
enormously increased, when we remember the extreme laboriousness with which 
the individual results would be obtained, when required, in the absence of a 
table. It would require, for example, a long day’s work to prove that 5,999,993 
is a prime, or even to resolve into factors 5,999,999, which has no factor less 
than 1000; in the latter case the verification would be simple enough, but in 
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the former the whole calculation would have to be duplicated before we could 
put any confidence in the result. In fact the contrast between the difficulty of 
obtaining the individual results without a table and the comparative ease with | 
which the table, as a whole, is constructed is most striking. | 
A portion of Burckhardt’s introduction relates to the calculation of the 
| 


logarithms of large numbers by means of a factor table in connection with such 
tables as Wolfram’s table of the hyperbolic logarithms of primes up to 10,009 
? to 48 decimal places; and the introduction to the fourth million contains an ex- | 


ample, the principle being to multiply by a factor whose logarithm is known, such 


q that the product shall be within the range of the factor table, and shall on trial 
be found to differ by only one or two units from a composite number, +, all of 
whose prime factors are within the range of the logarithmic table. The loga- ; i 
rithm is then obtained by the use of the formula : 


a® ° 
log (+ + a) =log x + -—3-—, —ete., 


in which, of course, the object is to render 7 as small as possible. *. 


‘* But,” to quote Mr. J. W. L. Glaisher (Camb. Phil. Proc., Feb. 11, 1878), “the 
real importance of an extended factor table lies in the fact that such a table is a 
fundamental one in the Theory of Numbers. The number of factors of a num- 
ber, their sum, etc., are elements which enter into many results in the Theory of 
Numbers, and it is clear that, even on this account alone, it is desirable to have 
a table by means of which the resolution of high numbers into their factors is 
rendered practicable, so that suitable examples, verifications, etc., of such theo- 
rems may be obtainable. Also, conjectural theorems of this nature could not 
readily be tested without such tables. Again, the law of frequency of prime 
numbers has been the subject of analytical investigation, and theoretic formula 
giving the average frequency of primes, and the approximate number of primes 


between any given limits, have been obtained by Legendre, Tchebycheff, Har- 
greave, Riemann, and others, and it is essential to have the means of comparing 


t the numbers given by these formula with those found by actually counting the 
primes.” 
AL Legendre published in 1808 his formula 


log + — 1.08366 
for the number of primes inferior to x. The reasoning by which he arrives at 


the form 
x 


A log « —B 
is vague and unsatisfactory, and no account is given of the method by which 
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the constants were determined. Legendre gives comparisons of the value of 
this formula with the number of primes counted, up to 400,000, from Vega’s 
Jabule. WLegendre’s confidence in this formula was perhaps based upon the 
fact that, although the constant may have been determined from the value 
= 10,000, the results were satisfactory up to 400,000. He says, ‘‘II est im- 
possible qn’ une formule représente plus fidélement une série de nombres d’ une 
aussi grande étendue et sujette nécessairement a de fréquentes anomalies.”’ 
Subsequently Legendre extended his comparisons to 1,000,000 by means of 
Chernac’s Crzbrum and the results were still satisfactory. 

In 1849 Gauss refers, in a letter to Encke, to the value of the constant in 
Legendre’s formula. He had, with the assistance of Goldschmidt, made an 
enumeration of primes from Burckhardt's three millions. He gives comparisons 
with Legendre’s formula, and calculates the values of the constant which must 
be used in order to bring the formula into agreement with the numbers counted 
for + == 500,000, T¥,000,000, etc., to 3,000,000. It appears that the constant must 
be decreased as yx increases. It was subsequently proved by Tchebycheff that 
the constant must be unity when + is infinite. Gauss also gave comparisons 
with the logarithmic integral 


ax 
li = 
9 log x 
Tchebycheff in 1848 arrived independently at the formula P 


dx 
log x 


and Hargreave in the following year represented the number of primes between 


x and by 


In 1872, Mr. J. W. L. Glaisher, in a paper read before the British Association, 
gave the result of enumerations made from Burckhardt’s and Dase’s Tables, and 
comparisons, throughout the second and ninth millions, with the differences of 
the values of the integral li +. The enumerations of Gauss were found to be 
very inaccurate. In 1876 and 1877, in communications to the Cambridge 
Philosophical Society, Mr. Glaisher, having carefully revised the enumerations, 
gives the results for the six millions which had then been published, and states 
his purpose of comparing the results with Legendre’s formula and with the 
logarithmic integral throughout the six millions. These enumerations were 
extended to the fourth and fifth millions, and the results appeared in the intro- 
ductions to those volumes; and the introduction to the sixth million contains tables 
giving the number of primes in each successive group of 1,000 numbers from 
0 to 9,000,000. Also in each group of 10,000 numbers, and in groups of 100,- 
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000, of 250,000, of 500,000, and in each million. The irregularity of the dis- 

tribution of primes is strikingly shown by the differences which are added in 

the tables where the groups contain 100,000 numbers or more; ‘‘not even in 

the table where the groups contain 500,000 numbers do the numbers of primes 

in the successive groups always decrease, there being an increase in one instance.” 
The results for the complete millions are as follows: 


Number of primes Number of primes 
# inthezth million. Difference. inferior to # millions. 
I 78,499 8066 78,499 
70,433 148,932 
3 67,885 216,817 
4 66, 329 960 283,146 
6 64, 336 —— 412,851 
7 63,799 hea 476,650 
8 63,158 68 539,808 
9 62,760 39 602,568 


In 1859 Riemann communicated to the Berlin Academy a paper (which 
appears at page 136 of his collected works) in which he investigates the law of 
frequency of primes, and finds for the non-periodic part of the expression for 
the number of primes inferior to x the formula 


re 
in which the terms are of the form + —- lirn, where # is the product of 


primes all different, and the sign is -+ or —, according as the number of primes 
is even or odd. To this memoir the late Professor H. J. S. Smith referred in 
1876 as containing ‘‘the only investigation of the asymptotic frequency of the 
primes which can be regarded as rigorous.”’ 

The introduction to the sixth million contains the complete comparisons, 
not only with Legendre’s formula, and with li +, but also with Riemann’s formula 
at intervals of 50,000 for the whole nine millions. The tables also give the 
deviations of the formulas from the actual numbers of primes. These deviations 
show the great superiority of Riemann’s formula, which is found to be more ac- 
curate than Legendre’s even for the smaller values of x, and to give satisfactory 
results for the whole nine millions. Legendre’s formula, on the other hand, 
commences at 650,000 to give results which are, with a single exception (at 
1,200,000) always in excess, the excess amounting to 391 at 8,450,000 and to 
282 at 9,000,000. The formula li x gives results largely in excess, even for the 
smaller values of x, and steadily diverges, but not so rapidly as Legendre’s 
formula; so that it is more accurate than the latter when x exceeds five millions. 
The mean deviations derived from the table are, for Riemann’s formula, — 9, 
for li x, + 163, and for Legendre’s formula, + 171. 
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To exhibit the deviations more clearly to the eye, Mr. Glaisher has given 
a diagram in which successive values of the mean of three deviations in the 
table are taken as ordinates corresponding to the middle values of # as abscissas. 
These ‘deviation curves” are, of course, extremely irregular, owing to the 
anomalous distribution of the primes; but, while the other curves tend steadily 
upward, Riemann’s curve shows no tendency to cumulative error; in fact it cuts 
the axis 18 times. Had the actual values of the deviations at intervals of 
50,000 been constructed as ordinates, the curve would have been still more 
irregular, but it would have cut or touched the axis 51 times within the nine 
millions, the last time at 8,500,000; the deviation at 9,000,000 is — 132. 

It is possible to determine the number of primes, for particular values of +, 


without actually constructing the factor table. The number of multiples of £, 


including / itself, inferior to + is F @ , the greatest integer which does not 


2 
exceed Ps Taking account of the numbers divisible by both f, and f,, the 


number divisible by either orf, is (=) + fa —E E | ; and in 
2 1 P2 


general the number divisible by one of the primes f,, fy, Py iS 
1 Pi Ps? 


The employment of this formula it will be noticed is equivalent to calculat- 
ing the number of entries of one of the primes f,, . . . f, that would be made 
by the ‘‘sieve”’ method, if a factor table were to be constructed to extend to +, 
but in the case of primes greater than f+ which occur as least factors of two- 
factor numbers only, within the limit of the table, the number of times p would 
be entered may be determined by a method analogous to the ‘‘ multiple’? method 


described above, 7@. ¢., it is only necessary to count the primes between / and rs 


Thus the whole number of entries may be determined, and by subtraction the 
number of primes. The process is, of course, very laborious, but Meissel in 
1870 and 1871 employed it to calculate the number of primes inferior to one 
million, to ten millions, and to one hundred millions. Meissel’s number for one 
million was exactly confirmed by the actual count from Burckhardt’s table, and 
Mr. Glaisher regards his results for ten millions and one hundred millions, 
namely, 664,580 and 5,761,461, as entitled to be accepted with confidence. 
Accordingly the values of the several formula were calculated for 2 — 10,000,000 
and for == 100,000,000. At ten millions, Riemann’s formula gives an ex- 
cess of 87, so that the deviation curve crosses the axis at least once between nine 
and ten millions; the deviation of li + is 338, and that of Legendre’s formula is 
560. At one hundred millions Riemann’s deviation is only go, that of li «a is 


} 
{ 
‘ 
i 
ne 
ps 
by 


EXERCISES. 23 


748, and that of Legendre’s formula is 6543. The accuracy of Riemann’s 
formula at this point is extraordinary, the error being less than gy },j)9 of the 
actual number of primes. 

The introduction to the sixth million, which contains many important and 

interesting results besides those briefly summarized in this review, concludes with 
a discussion of the mode’ of calculating the values of li + and of Riemann’s 
formula. 
, A mathematical table of this magnitude and importance may well be 
reckoned one of the ‘‘long results of time;’’ and England may be justly 
proud that such a work has not only been completed by one of her own mathe- 
maticians, but furnished with so worthy a crown. 


EXERCISES. 


1 
A HORIZONTAL WIND blows against a hemispherical dome of radius 2’. The 
pressure exerted by the wind on a plane surface normal to its direction is P 
pounds to the square foot; on a surface oblique to its direction the pressure ex- 
erted is normal, but is reduced in the ratio(Poncelet, Wecanique Industrielle, 403) 
where z is the angle of incidence. It is required to find the magnitudes and 
points of application of the horizontal and vertical components of the resultant 
wind-pressure. [ W. M. Thornton. | 
2 
THREE Closely connected tanks, 7,, 7,, 73, contain Q, gallons of water, Q, 
gallons of vinegar, Q, gallons of brandy, respectively. A flow is set up from 
7, through 7, to 7, and back to 7, at the rate of 1 gallon per second. The 
liquids are assumed to mix instantaneously, and the lengths of the connecting 
pipes are neglected. Show how to calculate the amount of water in each tank 


at the end of 7’ seconds. [ De Volson Wood. | 
3 
* To find the radius of a sphere, which being let fall into a given conical 
wine-glass, full of water, will displace the maximum or minimum quantity of 
water. [ Benjamin Alvord. | 


On one of the bounding radii of a given quadrant, radius 7, a semicircle is 
drawn, radius $7, the semicircle being within the quadrant. Find the average 
area of the circle touching the other bounding radius of the quadrant and the 
arc of the semicircle. [ Artemas Martin. | 
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5 
THREE LINEs are drawn from a point within a square to three corners of 
the square. Find a side of the square. [ J. M. Greenwood. | 
6 


F inp the value of the definite integral 


1 
(1 —} cos? #+-sin sin 27) 
[Z. G. Barbour. 
PUBLICATIONS RECEIVED. 


[In an early number a review will be given of the latest works on Algebra. ED.] 


A PracricaL AritHMetic. By G. A. Wentworth and Thomas Hill. Boston: Ginn, 
Heath, & Co. 

ELEMENTS OF ALGEBRA. By G. A. Wentworth. Boston: Ginn, Heath, & Co. 

ELEMENTS OF PLANE AND SoLip GeometTRY. By G. A. Wentworth. Boston: Ginn, 
Heath, & Co. 

PLANE AND SPHERICAL TRIGONOMETRY. By G. A. Wentworth. Boston: Ginn, 
Heath, & Co. 

SURVEYING AND Navication. By G. A. Wentworth. Boston: Ginn, Heath, & Co. 

A CoMPLETE ALGEBRA. By A. Schuyler. Cincinnati: Van Antwerp, Bragg, & Co. 

ELEMENTS OF GeoMETRY. By A. Schuyler. Cincinnati: Van Antwerp, Bragg, & Co. 

Five-PLACE LOGARITHMIC AND TRIGONOMETRIC TaBLes. Arranged by G. A. Wentworth 
and G. A. Hill. Boston: Ginn, Heath, & Co. ° 

A PRACTICAL ARITHMETIC. By Edward Olney. New York: Sheldon & Co. 

THe Compete ALGEBRA. By Edward Olney. New York: Sheldon & Co. 

ELEMENTARY GEOMETRY. New Edition. By Edward Olney. New York: Sheldon 
& Co. 

University ALGEBRA. By Webster Wells. Boston: Leach, Shewell, & Sanborn. 

ELEMENTS OF SURVEYING AND LEVELING. By Charles Davies. Revised by J. Howard 
Van Amringe. New York: A. S. Barnes & Co. 

THe THEORY oF Tursines. Abstract. By Robert H. Thurston. (Journal of the 
Franklin Institute. ) 

ALIGNMENT CURVES ON ANY SURFACE, WITH SPECIAL APPLICATION TO THE ELLIPSOID. 
Abstract. By C. H. Kummell. (Philosophical Society of Washington. Mathe- 
matical Section.) 

Tue THEORY OF ERRORS PRACTICALLY TESTED BY TARGET-SHOOTING. Abstract. By 
C. H. Kummell. (Philosophical Society of Washington. Mathematical Section. ) 

Batuistics. By James M. Ingalls. (Course of Artillery.) Fort Monroe: United 
States Artillery School. 

On CERTAIN PossiIBLE ABBREVIATIONS IN THE COMPUTATION OF THE LONG-PERIOD IN- 
EQUALITIES OF THE Moon’s MovioN DUE TO THE DIRECT ACTION OF THE PLANETS. 
By G. W. Hill. (American Journal of Mathematics.) 


CompouND DererRMINANTS. By C. A. Van Velzer. (American Journal of Mathemat- 
ics. ) 
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